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In spite of the relevance of the proposal introduced in the recent work A. Abedi, N. T. Maitra 
and E. K. U. Gross, J. Chem. Phys. 137, 22A530, 2012, there is an important ingredient which is 
missing. Namely, the proof that the norms of the electronic and nuclear wavefunctions which are 
the solutions to the nonlinear equations of motion are preserved by the evolution. To prove the 
conservation of these norms is precisely the objective of this Comment. 



In a remarkable recent work [1], Abedi et al. present 
an exact factorization of the molecular wavefunction into 
a nuclear and an electronic part, which allows to rig- 
orously introduce generalized and very useful concepts, 
such as the time-dependent potential energy surface. 
This formalism also sets the stage to better understand, 
and hence probably improve, very much used quantum- 
classical schemes, such as Ehrenfest, surface-hopping or 
Born-Oppenheimer dynamics. 

In spite of the relevance of the proposal, we consider 
that there is an important ingredient which is missing. 
Namely, the proof that the norms of the two functions 
$R(r, t) and x®!^)? which are solutions of the nonlin- 
ear Eqs. (28) and (29) in [1]. This is a key point in order 
to associate $R(r, t) and x(S:, t) to a marginal and a con- 
ditional probability amplitude, respectively, thus leading 
to their identification as nuclear and electronic wavefunc- 
tions, as it is the purpose of [1]. 

If one wants to be sure that a given set of equations of 
motion do conserve any given quantity, there are essen- 
tially two options: Either one explicitly forces the con- 
servation at the action level, e.g., using Lagrange multi- 
pliers, or one shows that there is another reason (e.g., a 
symmetry of the action) why the obtained equations of 
motion produce the conservation. In [1], neither of these 
two things are explicitly done. 

In these works it is proved that the exact solution, 
^'(R, r, t), of the time-dependent molecular Schrodinger 
equation can be written as a single product of the form 

*(E,E,0 =%fci)x(l,i) , (1) 
such that the partial normalization condition (PNC), 

y"dl|%(E,t)|2 = l VS,i, (2) 

is satisfied. This condition implies that also x is normal- 
ized if 4* is. 

But the PNC is only used in [1] to simplify the Euler- 
Lagrange equations once they have been obtained from 



the stationary action principle. As we said, since the 
conservation has not been enforced at the action level 
from the very first step, this occasional use of the PNC 
does not guarantee, in principle, that it holds for all times 
if no further proof is provided. 

But, before detailing our proof of this property does 
hold, let point out that two possible ways of proving it 
have been discarded here for different reasons. First, one 
could have shown that the action has a certain symme- 
try and obtain the conservation law as an application of 
Noether's theorem. We have been unable to find such a 
symmetry. Second, notice that the equations of motion 
in [1] can be written as 

idt^^(x, t) = X, 3t%] $g^(£, t) , (3a) 

idtx{R,t) - Hx[%X,5t%]x®,0 . (3b) 

If the operators and were linear and Hermitian, 
the conservation of the norm of the functions x 
would be straightforward. As that is not the case, a 
more careful analysis is in order. In the literature, some- 
times this type of problems are solved in a formal way. 
For some nonlinear cases, this formal procedure is useful 
and, assuming that all necessary conditions on the corre- 
sponding infinite dimensional Hilbert spaces are satisfied, 
it allows us to prove the conservation of the norm of the 
functions in a very simple way. However, in this particu- 
lar case, and due to the dependence of the time derivative 
of "I>R(r, t) in the definition of the operators, the only way 
to have Hermiticity is to show in advance that the norm 
of <I>R(r, i) is indeed conserved. This renders the reason- 
ing circular, and thus invalid. 

The problem of obtaining the dynamical equations can 
be addressed from a different perspective. Consider then 
a decomposition of the form given by Eq. (1) but without 
imposing the PNC condition. In this case, we define: 

f{R,t):^ [ di\^^i^,t)\' , (4) 
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and, if we use that is normalized at all times, we have: 
J dRf{R,t)\x{R,t)\^ = 1 . (5) 

If we consider this factorization for ^E" and introduce 
it in the variational framework used in [1], we obtain as 
dynamical equations: 



Ma 



+ ($R|Te + T„ +Ve + Wen " idtl^U.) 



X , 



(6) 



\x\Hf[Te+Tn+Ve + Wen] 



($r| (Te + T„ + K + Wen " i^t) |$r) - 



'i'R , (7) 



where the dependencies have been omitted, and 

Ve{l,t)--^Wee{r) + Ve%t{L,t) , 

K(E, t) : = WnniR) + KIM. t) , 

using the notation in [1]. 

These equations are completely general. Any factoriza- 
tion of the form (1), where the molecular wave function 
^ is a solution of the molecular Schrodinger equation, 
satisfies them. Also, notice that, if we have f(R,t) — 1 
for all t (i.e., the factorization satisfies the PNC for all 
values of time), then these equations reduce to Eqs. (28) 
and (29) in [1]. Now, the factorization in (1) exhibits 
an invariance under the group of invertible functions on 
Co = C - {0} given by 



a(i,i)x(l,i) , 



(8) 



a(R,t) 



where a(R,, t) is any complex function without zeros. In 
other words, for any given solution, x(Rj t) and $R(r, t), 
of the equations of motion in (6) and (7), we can obtain 
new solutions, xCRit) a-nd $R(r, i), which produce the 
same molecular wavefunction ^'(R, r, t), by applying the 
above transformation. Of course, these new functions 
will be solution to the equations with the corresponding 
/(R, t). Also notice that this gauge freedom enlarges the 
C/(l)-freedom discussed in [1], where only the phase of 
each function is transformed. 



Let us consider now a gauge fixing defined as 
a(R,t) = e^^(l'*\//(R,t) 



(9) 



where /(R, t) is defined by Equation (4) and 6'(R, t) is 
arbitrary. Notice that we are extending the usual notion 
of gauge fixing in a subtle way. The transformation is 
considered for the full trajectory ^'(R, r, i), since it de- 
pends explicitly on the norm of the function $r along 
it. Thus it becomes a dynamical property. In particu- 
lar, if we consider a solution with initial unit norm, i.e., 
/(B;7 0) = 1, VR, we find that the initial conditions for 
the original and the transformed curves coincide: 

x(l, 0) = x(E, 0) , %(£, 0) = %(i, 0) , (10) 

and that functions <f>R and x normalized by construc- 
tion for all values of time: 



dR|x(R,OI' = 1 = / dr||.R(r,0|' yt,R . (11) 



Finally, consider any solution of Eqs. (28) and (29) 
in [1] for some initial conditions x^lS) and 'I'Rd) that 

satisfy the PNC, and let us ask whether or not the 
PNC is satisfied at subsequent times. We have seen 
that, among the factorizations of the molecular wave- 
function 4'(R, r, i) [with initial conditions ^(R, r, 0) = 

X°(R)$R(r))] there is one, given by x(R, t) and l'R(r, t), 

which satisfies the PNC for all values of time, and also 
Eqs. (6-7) with /(R, t) — 1 and initial conditions x^lS) 
and ^'R(r). Now, given that Eqs. (6-7) with /(R,t) = 1 
are precisely Eqs. (28) and (29) in [1], wc have that the 
functions xCR^t) and <I>R(r, t) are a also a solution to 
them, with initial conditions x°(R:) and $^(r), and such 
that their norms are preserved for all time. Since we 
have not modified the initial conditions, if we assume 
that Eqs. (28) and (29) in [1] have unique solution for 
given initial conditions, we can conclude that the norm- 
conserving solution to them that we have found must be 
the arbitrary one we began with, thus showing that every 
solution of Eqs. (28) and (29) in [1] with initial conditions 
that satisfy the PNC also satisfies it at all times. 
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